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The Bethe-Ansatz local density approximation (LDA) to lattice density functional theory (LDFT) 
for the one-dimensional repulsive Hubbard model is extended to current-LDFT (CLDFT). The trans- 
port properties of mesoscopic Hubbard rings threaded by a magnetic flux are then systematically 
investigated by this scheme. In particular we present calculations of ground state energies, persistent 
currents and Drude weights for both a repulsive homogeneous and a single impurity Hubbard model. 
Our results for the ground state energies in the metallic phase compares favorably well with those 
obtained with numerically accurate many-body techniques. Also the dependence of the persistent 
currents on the Coulomb and the impurity interaction strength, and on the ring size are all well 
captured by LDA-CLDFT. Our study demonstrates that CLDFT is a powerful tool for studying 
one-dimensional correlated electron systems with high accuracy and low computational costs. 
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I. INTRODUCTION 

Quantum dots, routinely made by electrostatically 
confining a two-dimensional electron gas [T], have been 
extensively studied in recent years [2]. The interest in 
these low-dimensional structures stems from the fact that 
their physics is controlled by quantum effects. Further- 
more, while sharing many similarities with real atoms, 
quantum dots manifest intriguing low-energy quantum 
phenomena, which are specific to them. This is because 
their properties can be influenced by external factors such 
as the geometry or the shape of the confining potential 
and the application of external fields. Clearly some of 
these features are not accessible in real atoms. Research 
in the past has been motivated by the possibility of de- 
veloping novel quantum dot based devices in both the 
fields of quantum cryptography/computing [3] and spin- 
tronic [I], as well as by the simple curiosity of exploring 
the properties of many-electron systems in reduced di- 
mensions. 

Quantum rings represent a particular class of quan- 
tum dots [5,6, where electrons are confined in circular 
regions [7J |S] ■ The circular geometry can sustain an elec- 
trical current, which in turns can be induced by threading 
a magnetic flux across the ring itself. Such a magnetic 
flux produces exciting effects like Ahanorov-Bohm (AB) 
oscillations [DJ HPJ and persistent currents [TT], effects 
that were anticipated as early as the late 60's [T51 IT5] . 
In one dimension (ID) the persistent currents have been 
thoroughly studied [UltS]. These, as many other phys- 
ical properties of the ring, are a periodic function of the 
magnetic flux quantum, $o = hc/e (h is the Planck's con- 
stant, c the speed of light and e the electron's charge). 

A number of earlier theoretical studies [15TtT9] on per- 
sistent currents focused on unveiling the role of elec- 
tron correlations and disorder over the electron trans- 
port. This line of research is inspired by the fact that 
electronic correlation in ID always leads to non-fermionic 
low-energy quasiparticle excitations. In fact, even in the 



presence of weak interaction, ID fermions behave dif- 
ferently from a Fermi liquid and their ground state is 
generally referred to as Luttinger liquid. This possesses 
specific collective excitations [2D] . 

There are two theoretical frameworks commonly used 
to study finite ID rings [5T]. The first is based on the con- 
tinuum model, where electrons move in a uniform neu- 
tralizing positive background and interact via Coulomb 
repulsion, e 2 /Aire^r (eo is the vacuum permittivity, e the 
electron charge, r the distance between two electrons). 
The second is populated by lattice models, where the 
electronic structure is written in a tight-binding form and 
the electron-electron interaction is commonly described 
at the level of Hubbard Hamiltonian [55] [53]. In both 
frameworks exact diagonalization (ED) has been the pref- 
erential solving strategy for small systems (small number 
of sites and electrons) [HI E3] . Additional methods used 
to study quantum rings over lattice models include Bethe 
Ansatz (BA) [25, 26 , renormalization group [57] and den- 
sity matrix renormalization group |28j . In contrast the 
continuum model has been tackled with self-consistent 
Hartree Fock techniques [5DJ, bosonization schemes [3D], 
conformal field theory [31] . current-spin density func- 
tional theory [35] and quantum Monte Carlo [33 . 

Many of the methods developed for solving lattice 
models for interacting electrons suffer from a number of 
intrinsic limitations connected either to their large com- 
putational overheads or to the need of using a drasti- 
cally contracted Hilbert space. Density functional the- 
ory (DFT) can be a natural solution to these limitations. 
DFT is a highly efficient and precisely formulated method 
[331 ES] , originally developed for the Coulomb interaction 
(this is commonly known as ab initio DFT), and then ex- 
tended to lattice models [36ti38] , Lattice DFT (LDFT) is 
based on the rigorously proved statement that the ground 
state of an interacting electron system is a universal func- 
tional of the local site occupation. The functional, as in 
ab initio DFT, is unknown explicitly. However all the 
many-body contributions to the total energy can be in- 
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corporated in a single term, the exchange and correlation 
(XC) energy, for which a hierarchy of approximations can 
be constructed. 

The most commonly used approximation for the XC 
energy in ab initio DFT is probably the local density ap- 
proximation (LDA) [35| I39j , where the exact (unknown) 
XC energy is replaced by that of the homogeneous elec- 
tron gas. The theory is then expected to work best in 
situations close to those described by the reference sys- 
tem, i.e. close to the homogeneous electron gas. Since in 
ID Fermi liquid theory breaks down, the homogeneous 
electron gas is no longer a good reference. For the ho- 
mogeneous Hubbard model it was then proposed [35] HU] 
to use instead the BA construction of Lieb and Wu [26] . 
Such a scheme was then applied successfully to a wide 
range of situations |4"TH4"T] and more recently it has been 
extended to time-dependent problems [HO and to 

the 3D Hubbard model [5T] . 

LDFT can be further extended to include the action of 
a vector potential, i.e. it can be used to tackle problems 
where a magnetic flux is relevant. This effectively corre- 
sponds to the construction of current-LDFT (CLDFT). 
Such an extension of LDFT was proposed recently for 
one-dimensional spinless fermions with nearest-neighbor 
interaction [52j and it is here adapted to the repulsive 
Hubbard model. The newly constructed functional is 
then used to investigate total energies, persistent cur- 
rents and Drude weights of a mesoscopic repulsive Hub- 
bard ring threaded by a magnetic flux. 

The paper is organized as follows. Section II reviews 
the theoretical foundations leading to the construction 
of CLDFT and to its LDA. Then we present our results 
for both homogeneous and defective rings, highlighting 
the main capabilities and limitations of our scheme, and 
finally we conclude. 



II. THEORETICAL FORMULATION OF 
CURRENT LATTICE DFT 

Current-DFT (CDFT) is a generalization of time- 
dependent density functional theory [53] to include in 
the Hamiltonian an external vector potential [54] . In 
this case the theory is constructed over two fundamen- 
tal quantities, namely the electron density, n, and the 
paramagnetic current density, j p . The Hohenberg-Kohn 
theorem [31] is thus expanded to the statement that the 
ground state n and j p uniquely determine the ground- 
state wave-function and consequently the expectation 
values of all the operators [55] [55]. Equally important 
is the fact that the standard Kohn-Sham construction 
can also be employed for CDFT, so that the many-body 
problem can be mapped onto a fictitious single-particle 
one, with the two sharing the same ground state n and 
j p [55] [56]. Practically one then needs to solve self- 
consistently a system of single-particle equations. Also 
for CDFT all the unknown of the theory are incorporated 
in the XC energy, which then needs to be approximated. 



The scope of this section is to describe how ab initio 
CDFT has been translated to lattice models and how 
a suitable approximation for the XC energy associated 
to the Hubbard Hamiltonian can be constructed. Our 
description follows closely the one previously given by 
Dzierzawa et al. [52] . In general a vector potential, A, en- 
ters into a lattice model via Peierls substitution [57] [55] , 
where the matrix elements of the A-dependent Hamilto- 
nian, H(r,p+ -A), can be written in terms of those for 
A = as 

(R'\H{r,p+-A)\R) = (R'\H(r,p)\R)e-^fi (l) 
c 

where c is the speed of light and \R) is the generic orbital 
located at the position R and belonging to the basis set 
(here assumed orthogonal) used to construct the tight- 
binding Hamiltonian. 

When the Peierls substitution is applied to the con- 
struction leading to the ID Hubbard model the only term 
in the Hamiltonian that gets modified is the kinetic en- 
ergy T . This takes the form 

L 

f=-t J2^ ,/L ^li + i^i+hc), (2) 

ct, 1=1 

where we have considered a system comprising L atomic 
sites (note that the ring boundary conditions imply 
L + 1 = 1). In the equation |i| (c a i) is the cre- 
ation (annihilation) operator for an electron of spin a 
(cr =t, I) at the /-site, t is the hopping integral and Q a i 
is the phase associated to the l-th bond, which effec- 
tively describes the action of A. The remaining terms in 
the Hamiltonian are unchanged so that the ID Hubbard 
model in the presence of a vector potential is defined by 

L 

Hubbard = f + f/ + 5>f rt n n ( 3 ) 
I 

where {f ; cxt } is the external potential (vf Kt is the on- 
site energy of the Z-site), while the Coulomb repulsion 
term is U = U fi^ih^i , with U being the Coulomb 
repulsion energy and n a i = $ al c a i- Throughout this work 
we always consider the diamagnetic (non-spin polarized) 
case so that ^-j-/ = = and = = m. 

The first step in the construction of a CLDFT is the 
formulation of the problem in a functional form. The 
basic variables of the theory are the site occupation n\ — 

|tt.; and bond paramagnetic current, ji — (^\ji\^f), 
where |^/) is the many-body wavefunction and the bond 
paramagnetic current operator is defined as 

3l = -it(e- i *'/ L ct, l+1 c al -hc) (4) 

In complete analogy to ab initio CDFT we can write the 
total energy, £, of the Hamiltonian ^ as a functional of 
the local external potentials and phases 

£ = F[nuji] + X>r% + E ( 5 ) 
i i 
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so that 



m = (ni) = 



3 1 = Ol) = 



d£ 

9wj oxt 
dS 



(G) 



J-[ni,ji] is a universal functional, whose functional 
derivatives with respect to {n{\ and {ji} satisfy the fol- 
lowing two equations 



v, = - 



dni 
dJF_ 
dji ' 



(7) 



Note that equations ^ through Q follow directly from 
the properties of the Legendre transformation. 

In order to make the theory practical one has now to 
introduce the auxiliary single-particle Kohn-Sham sys- 
tem. This is described by a single-particle Hamiltonian, 
H s , whose ground state site occupations and bond para- 
magnetic currents are identical to those of the interacting 
system [described by equation (13b] . H s reads 



i 



VlUl, 



(8) 



where f s = -tJ2tj=i( e ~ l *' /L + he) and the 

associated local effective potentials and phases are vf and 
respectively. The single-particle Schrodinger equation 
is then 



and the site occupation is defined as 

n s t =J2UK\ni\K), 



(9) 



(10) 



where f a is the occupation number. An analogous ex- 
pression can be written for j ; s . 

The energy functional associated the Kohn-Sham sys- 
tem, J-" s , can be constructed by performing again a Leg- 
endre transformation 



(ii) 



where £ s is the total energy of the single-particle system 
and the following two equations are valid 



$! = - 



d 3l 



(12) 



The crucial point is that in the ground state the real and 
the Kohn-Sham systems share the same site occupation 
and paramagnetic current, i.e. ni = nf and ji = jf. 



Thus one is now in the position of defining the XC 
energy, £ xc , as usual, i.e. as the difference between T 
for the interacting and the Kohn-Sham systems after the 
classical Hartree energy £ H has also been subtracted, 



£ xc [n/, ji] = F[n u j{\ - J^[ni,ji] - £ H [ni] 



(13) 



Note that for all the functionals in equation ( 13 1 we took 
the short notation {n{\ — > ni and {ji} — ¥ ji, i.e. the func- 
tionals depend on all the on-site occupations and bond 
paramagnetic currents. The single-particle effective po- 
tentials and phases can now be defined. In fact by taking 
the functional derivative of equation ( |13[) w ith respect to 
ni and ji and by using the equations (17j) and ( 12 ) one 
obtains 



u l 



where 



(14) 



dm 



nexe 

1 dji ' 



(15) 



and vj 1 = dSi/dni (= Uni/2) is the Hartree potential. 

Finally £ xc can be re-written in terms of the expecta- 
tion values of the original Hamiltonian. In fact by substi- 
tuting the functional forms of JF and JF S into the equation 
PI 



I I 

(16) 

by using the equations ^ and 
$3(«f " wf**)"! = £ s - £ - (V B \f s |* s 



(&\t + u\v) , 

(17) 

and again by substituting equation ( fl7| ) into equation 
(16 1, one obtains a close expression for the XC energy 



g*c = (*|T + [/|*)-(* s |f s |* s )+£($^$ i )ji-^ H [^]- 

(18) 

Once the theory is formally established the remaining 
task is that of finding an appropriate approximation for 
£ xc . As for the case of standard LDFT gDJ 01], the 
strategy here is that of considering the BA solution for 
the homogeneous limit of ^Hubbard (^ ms ^ s defined in 
equation ^ by setting vi = v and = $) and then of 
taking its local density approximation n — > ni, $ — > 



i.e. 



1 



(19) 



where e xc [n,j) = £ KC [n, j]/ L is the XC energy density 
(per site) of the homogeneous system. The first term of 
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the equation ( 18 1 can be calculated exactly using the BA 
procedure |59j. This provides the ground state energy as 
a function of n and <&, so that one still needs to re-express 
it in terms of n and j. However the phase variable $ can 
be eliminated from the ground state energy by using 
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d£{n, $) 
9$ 



(20) 



Thus finally one can explicitly write e xc (n, j) (the full 
derivation for the ID Hubbard Hamiltonian is presented 
in the Appendix) 



e x >,j) = e x >,0)+ *A xc (n)j 2 , 



where 



e xc (n,0) 
A xc (n) = 



£ BA (n,0)-£°(n,0)-£ ll (n) 



L 



D°(n) £>ba(„) 



(21) 



(22) 



In the equations above £°(n,0) and D®(ri) are respec- 
tively the non-interacting ground state energy and charge 
stiffness, while £ BA (n, 0) and D BA (n) are the same quan- 
tities for the interacting case as calculated from the BA. 
Finally, the XC contributions to the Kohn-Sham poten- 
tial can be obtained by simple functional derivative (in 
this case by simple derivative) of the exchange and cor- 
relation energy density with respect to the fundamental 
variables n and j, i.e. 



7 BALDA 



{null) = 



de 



lc (n,j) 
dn 



(23) 



and 



BALDA 



(nidi) 



de xc (n,j) 
~~dj 



-A xc (n)j 



>rii,j->ji 



(24) 

stands for Bethe Ansatz local 



where BALDA, as usual, 
density approximation. 

In the two panels of figure [l] we present e xc (n,0) and 
A xc (n) as a function of the electron density, n, for dif- 
ferent interaction strengths U/t. As in the case of stan- 
dard LDFT also for CLDFT there is a divergence in the 
n-derivative of both e xc (n, 0) and A xc (n) at half-filing 
(n = 1). This is in correspondence of the metal-insulator- 
transition present in the ID Hubbard model for finite 
U jt. In the case of A xc (n) the divergence is also in A xc (n) 
itself. 

The solution of the Kohn-Sham problem proceeds as 
follows. First an initial guess for the site occupations 
is used to construct the initial local paramagnetic cur- 
rent density. Then, the functional derivatives of equa- 
tions (23) and (24) are evaluated at these given n and 
j so that the Kohn-Sham potential is constructed. The 
Kohn-Sham equations are then solved to obtain the new 
set of Kohn-Sham orbitals from which the new orbital oc- 
cupations and bond paramagnetic currents are calculated 




FIG. 1: (Color online) The XC energy density (per site) for a 
homogeneous ID Hubbard ring threaded by a magnetic flux 
as a function of the electron density and for different values 
of interaction strength U/t: (a) e xc (n,Q) and (b) A xc (ri). 



[by using equation ( |10| )]. The procedure is then repeated 
untill self-consistency is reached, i.e. until the potentials 
(or the densities) at two consecutive iterations vary be- 
low a certain threshold. After convergence is achieved 
the total energy for the interacting system is calculated 
from 



where the first term is the sum of single-particle energies 
and the other terms are the so-called double counting 
corrections. 



III. RESULTS AND DISCUSSION 

We now discuss how CLDFT performs in describing 
both the energetics and the transport properties of ID 
Hubbard rings in presence of a magnetic flux. For small 
rings our results will be compared with those obtained 
by diagonalizing exactly the Hamiltonian of equation (|3| , 
while CLDFT for large rings will be compared with the 
BA solution. First we will consider homogeneous rings 
and then we will explore the single impurity problem. 



A. Homogeneous rings: general properties 

In this section we focus our attention on discussing 
the general features of CLDFT applied to homogeneous 
Hubbard rings threatened by a magnetic flux, i.e. on 
the performance of CLDFT in describing the Ahanorov- 
Bohm effect. We start our analysis by comparing the 
CLDFT results with those obtained by ED. Since ED 
is numerically intensive such a comparison is limited to 
small systems. 

In figure [2] we present the first low-lying energy levels, 
£ , calculated by ED as a function of the magnetic flux, 
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<f>, for a small 12-site ring at quarter filling (n = 1/2). In 
particular we present results for the non-interacting case 
[panel (a)] and for the interacting one at three different 
interaction strengths: (b) U/t — 2, (c) U/t — 4 and (d) 
U/t = 6. Exact results (ED) are in black, while those 
obtained with CLDFT in red. In general the ground 
state energy is mimimized at $ = when the number 
of electrons is N = 4m + 2 and at $ = ir for N = 4m, 
with m being an integer [60j . Here we consider the case 
TV = 4m + 2 where the ground state is a singlet [61] . 

For non-interacting electrons, U/t = 0, the total en- 
ergy of the singlet ground state is a parabolic function of 
Also the various excited states have a parabolic de- 
pendence on $ and simply correspond to single-particle 
levels with different wave-vectors. As the electron- 




-2-10 1 2-2-10 12 

<t>/7t <E>/ji 



point that the CLDFT energy is practically identical to 
that calculated with ED. However CLDFT completely 
misses the cusps in the £ ($) profile arising from the 
crossover between the singlet and the triplet state. Level 
crossing invalidates the BA approximation leading to the 
interacting XC energy [see equation p6| ) in the appendix] 
and so failures are expected [52;. This observation is in 
agreement with earlier studies [21j in which the inability 
of CDFT to reproduce level crossing was already noted. 
Nevertheless, as long as the singlet remains the ground 
state, the agreement between CLDFT and ED results is 
remarkable, even if this small ring is rather far from be- 
ing a good approximation of the thermodynamic limit 
(the BA solution) upon which the functional has been 
constructed. 
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FIG. 2: (Color online) The low lying energy spectrum, £, of 
a 12-site ring at quarter filling (n — 1/2) as a function of 
the magnetic flux, $, and calculated for different interaction 
strengths U /t. The black dotted lines represent ED results 
while the dashed red ones are for CLDFT. Note that for the 
non interacting case, U/t = 0, in panel (a) there is no differ- 
ence between CLDFT and ED. Panels (b)-(d) are for the in- 
teracting case at different interaction strengths: (b) U/t — 2, 
(c) U/t = 4 and (d) U/t — 6. In panel (b) the blue arrows in- 
dicate the region where the triplet state becomes the ground 
state. 

electron interaction is turned on the non-interacting spec- 
trum gets modified in two ways. Firstly there is a second 
branch in the ground state energy as a function of $ ap- 
pearing at around $ = ±7r (see the blue arrows in panel 
(b) of Fig. [2]) . This originates from the degeneracy lifting 
between the single and the triplet solution at $ = ±7r, 
with the triplet being pushed down in energy and becom- 
ing the ground state. The $ region where the ground 
state is a triplet widens as the interaction strengths in- 
creases. The second effect is the expected reduction of 
the ground state total energy as a function of U/t. 

Since CLDFT is a ground state theory, it provides ac- 
cess only to the ground state energy, £ . This is calculated 
next and plotted in figure [2] in the interval — n < <I> < 7r 
for different U/t. As one can clearly see from the figure 
the performance of CLDFT is rather remarkable, to a 



FIG. 3: (Color online) Persistent current profile, j, for a 12- 
site ring at quarter filling (n = 1/2) obtained with both ED 
and CLDFT for different U/t. The full lines are the j calcu- 
lated with ED while the dashed ones are for CLDFT. 

Having calculated the total energies with both ED and 
CLDFT, the corresponding persistent currents, j, can 
be obtained by taking the numerical derivative of £ ($) 
with respect to $. In figure [3] we show results for the 
12 site ring at quarter filling (n = 1/2), whose total en- 
ergy was presented in figure [2j In particular we plot j 
only over the period — it < <f> < w, since all the quanti- 
ties are 2-7T periodic. The figure confirms the linearity of 
the persistent currents with the magnetic flux for all the 
interaction strengths considered. The same is also true 
for other fillings for the 12 site ring (not presented here) 
away from half-filling. We also observe that the magni- 
tude of persistent currents reduces with increasing U/t 
for both ED and CLDFT and that the precise depen- 
dence of j on U /t is different for different fillings. This 
is in good agreement with previous calculations based on 
the BA technique [63 . 

ED is computationally demanding and cannot be per- 
formed beyond a certain system size. For this reason, 
in order to benchmark CLDFT for larger rings, we have 
calculated the ground state energy with the BA method. 
An example of these calculations is presented in figure 
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|1J where once again we show £($) for L — 20, U/t = 4 
and different numbers of electrons. Also in this case the 
agreement between the BA results and those obtained 
with CLDFT is remarkably good as long as the ground 
state is a singlet. Interestingly we note that the agree- 
ment is better for low filling but it deteriorates as one 
approaches the half- filling case (TV — 20 in this case). 
This is somehow expected given the discontinuity of A xc 
and of the derivative of e xc at n — 1 (see figure [TJ, leading 
to the Mott transition. The presence of these disconti- 
nuities, although qualitatively correct, poses numerical 
problems and losses in accuracy. 
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FIG. 4: (Color online) Ground state energy, £($), as a func- 
tion of the magnetic flux, <E>, calculated with both the BA 
technique (dotted black line) and CLDFT (dashed red line). 
Calculations are carried out for L = 20, U/t — 4 and differ- 
ent numbers of electrons: (a) TV = 2 (n = 1/10), (b) TV = 6 
(n = 0.3), (c) TV = 10 (n = 1/2) and (d) TV = 14 (n = 0.7). 
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FIG. 5: (Color online) Drude coefficient D c as a function of 
the interaction strength U/t (top panel) and of the number of 
sites in the ring, L (bottom panel). All the calculations are 
for quarter filling and the results in the top panel are for a 
60-site ring. In the figure we compare CLDFT results (dotted 
black lines) with those obtained by the BA technique in the 
thermodynamic limit (dashed red lines). Calculations in the 
lower panel are for U/t = 2. 



The final quantity we wish to consider is the charge 
stiffness or Drude weight, D c , defined as 



D c = 



L d 2 £(n,< 5>) t 

2 



$=0 



(26) 



This is essentially the slope of the persistent current as 
a function of $ calculated at $ = and defines the mag- 
nitude of the real part of the optical conductivity in the 
long wave-length limit (see appendix for more details). 
D c determines both qualitatively and quantitatively the 
transport properties of the ring. Importantly in the limit 
of large rings it exponentially vanishes for insulators, 
while it saturates to a finite value for metals. Many stud- 
ies have been devolved to calculating D c for interacting 
systems. Romer and Punnoose have studied D c for fi- 
nite Hubbard rings using an iterative BA technique |62j . 
Eckern et. al. explored the relation between D c and the 
so-called phase sensitivity, A£, for spinless fermions. A£ 
is the difference in the total energy calculated at $ = 
(periodic ground state) and that at $ = tt (antiperiodic 
ground state) [64, 65j. Recently a density matrix renor- 
malization group algorithm has been developed to deal 
with complex Hamiltonian matrices and used to calculate 
D c for spinless fermions [28] , 



Since the agreement between CLDFT and ED is proved 
for small rings (the slopes of the persistent currents as a 
function of calculated with CLDFT and ED are essen- 
tially identical in figure[3| we concentrate here on a larger 
system, namely a homogeneous 60 site ring at quarter 
filling. Our results for the Drude weight as a function of 
U/t are presented in figure 5} Again the CLDFT data 
are compared with those calculated with the BA in the 
thermodynamic limit (L — > oo) and the agreement is 
rather satisfactory. We note that, as for the ground state 
energy, also for the Drude weight the CLDFT seems to 
perform less well as U/t increases, i.e. as the interaction 
strength becomes large. Then in the lower panel of figure 
[5] we illustrate the scaling properties of D c as a function 
of the number of sites in the ring, L (we consider quarter 
filling and U/t = 2). Clearly D c does not vanish at any 
lengths demonstrating that the system remains metallic. 
Furthermore it approaches a constant value already for 
L > 40. In the picture we also report the asymptotic 
value predicted by the BA in the thermodynamic limit 
L — > oo for this set of parameters. We find that the cal- 
culated CLDFT value is only 0.06% larger than the BA 
one, i.e. it is in quite remarkable good agreement. 
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B. Scaling properties 

Next we take a more careful look at the scaling prop- 
erties of the persistent currents and the Drude weights 
as a function of both the ring size and the interaction 
strength. It is well known that j is strongly size depen- 
dent, since it originates from electron coherence across 
the entire ring [55]. For a perfect metal one expect j to 
scale as 1/L [21] . In Figure [6] the value of the persis- 
tent currents as a function of the ring size are presented 
for different electron fillings and for the two representa- 
tive interaction strengths of U/t = 2 (a) and U/t = 4 
(b) . Calculations are performed with both the exact BA 
and CLDFT. As a matter of convention we calculate the 
persistent currents at $ = tt/2. 

In general we find a monotonic reduction of the persis- 
tent current with L and an overall excellent agreement 
between the BA and the CLDFT results over the entire 
range of lengths, occupations and interaction strengths 
investigated. A non- linear fit of all the curves of figure [6] 
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FIG. 6: (Color online) Persistent current, j, as a function of 
the number of site in the ring, L, and for different electron 
occupations, n: (a) U/t = 2, (b) U/t = 4. Results are ob- 
tained with both the exact BA technique and CLDFT. In the 
figure the persistent currents are calculated at $ = 7i"/2, i.e. 

i=j(T/2) 

returns us an almost perfect 1/L dependence of j with 
no appreciable deviations at any n or U /t. This indi- 
cates a full metallic response of the rings in the region of 
parameters investigated, thus confirming previous results 
obtained with the BA approach [63] . 

Then we look at the dependance of j and D c on the 
interaction strength. In this case we consider a 60 site 
ring and three different different electron fillings. In gen- 
eral for small fluxes one expects j = 2D C 4> and our nu- 
merical results of Fig. [7] demonstrates that this is ap- 
proximately correct also for our definition of persistent 
currents [j = j($ = 7r/2)] over the entire U/t range in- 
vestigated. We find that both j and D c monotonically 
decrease as a function of the interaction strength, essen- 
tially meaning that the predicted long-wavelength opti- 
cal conductivity is reduced as the electron repulsion gets 
larger. 

Also in this case the agreement between the BA and 
the CLDFT results is substantially good, although sig- 
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FIG. 7: (Color online) Persistent current, j, and Drude 
weight, D c as a function of interaction strength U/t for a 
60 site ring at different fillings. Results are obtained with 
both the exact BA technique and CLDFT. 



nificant deviations appear in the limit of large U /t and 
electron filling approaching half-filling. This again corre- 
sponds to a region of the parameter space where the XC 
potential approaches the derivative discontinuity. 

It was numerically demonstrated in the past |63| that 
the persistent current (and so the Drude weight) at half- 
filling follows the scaling relation j ~ c^ 11 '^, with £ ~ 1. 
However, to the best of our knowledge, no scaling relation 
was ever provided in the metallic case. We have then 
carried out a fitting analysis (the fit is limited to values 
of j and D c for U/t > 2) and found that our data can be 
well represented by the scaling lows 



j = jo{u/t)-e, 



D r 



D (U/t)-i 



(27) 



In general and as expected we find /? = 7 and a quite 
significant dependence of the exponents on the filling. 
In particular table [I] summarizes our results and demon- 
strates that the decay rate of both the persistent currents 
and the Drude weights increases as the filling approaches 
half-filling. Furthermore the table also quantifies the 
differences between the BA and the CLDFT solutions, 
whose exponents increasingly differ from each other as 
the electron filling gets closer to n = 1 (for n = 0.7 we 
find /3 BA ~ 2/3 CLDFT ). 



n 






0.3 


0.036 


0.036 


0.5 


0.085 


0.104 


0.7 


0.151 


0.246 



TABLE I: Exponents for the empirical scaling laws of equa- 
tion (271 as fitted from the data of figure [7] 



Finally, by combining all the results of this section we 
can propose a scaling law for both the persistent currents 
and the Drude weights, valid in the metallic limit of the 



Hubbard model, i.e. away from half-filling. This reads 



join) (U 



L 



-000 



(28) 



where both the constant j and the exponent (3 are func- 
tion of the electron filling n. Note that an identical equa- 
tion holds for D c . 



C. Scattering to a single impurity 

Having established the success of the BALDA to 
CLDFT for the homogeneous case we now move to a 
more stringent test for the theory, namely the case of 
a ring penetrated by a magnetic flux in the presence of 
a single impurity. Such a problem has already received 
considerable attention in the past [321 HS1 EZ] ■ Note that, 
as in ab initio DFT, this is a situation different from the 
reference system used to construct the BALDA (since it 
deals with a non homogeneous system) and therefore one 
might expect a more pronounced disagreement with ex- 
act results. As the BA equations are integrable only for 
the homogeneous case we now benchmark our CLDFT 
results with those obtained by ED. This however limits 
our analysis to small rings. 

The single impurity in the ring is described by simply 
adding to the Hamiltonian of equation ([3| the term 



(29) 



where £j mp is the modification to the on-site energy at 
the impurity site i. The inclusion of an impurity pro- 
duces in general electron backscattering so that we ex- 
pect the persistent currents to get reduced. In figure [8] 
we present the general transport features for this inho- 
mogeneous system. Calculations have been carried out 
with CLDFT for a ring comprising 53 sites and TV = 26, 
U/t = 4. Again the persistent currents are calculated at 
$ = tt/2. 
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Panel (a) shows j as a function of the impurity on- 
site energy. As expected from standard scattering theory 
the current is reduced as £; mp increases, thus creating 
a potential barrier. The electron density profile for this 
situation is presented in panel (b), where one can clearly 
observe an electron depletion at the impurity site and 
Friedel's oscillations around it. 

A quantitative assessment of our CLDFT results is pro- 
vided in Fig. [9] where they are compared with those ob- 
tained by exact diagonalization for a 13 site ring close 
to quarter filling (TV = 6). In particular we present j 
as a function of the impurity potential, £i mp , for both 
U/t = 2 and U/t — 4. In general we find a rather satis- 
factory agreement between CLDFT and the exact results 
in particular for small £i mp and U /t. As the electron scat- 
tering becomes more significant deviations appear and 
the quantitative agreement is less good. Importantly we 
notice that the ED results systematically provide a per- 
sistent current lower than that calculated with CLDFT, 
at least for the values of electron filling investigated here. 
This seems to be a consistent trend also present for the 
homogeneous case (see figure [7]), although the deviations 
in that case are less pronounced (for the same electron 
filling and interaction strength) . Therefore we tentatively 
conclude that most of the errors in the impurity problem 
have to be attributed to the errors already present in the 
homogeneous case. We then expect that CLDFT pro- 
vides a good platform for investigating scattering prob- 
lems at only minor computational costs. As such CLDFT 
appears as the ideal tool for investigating the interplay 
between electron-electron interaction and disorder in low 
dimensional structures. 
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FIG. 9: (Color on line) Comparison between the persistent 
currents calculated with CLDFT (black symbols and dotted 
line) and by ED (red symbols and dashed line) for a 13 site 
ring and TV = 6. The j's are obtained at $ = tt/2 for two 
different values of the interaction strength, namely U/t = 2 
(a) and U/t = 4 (b). 



FIG. 8: (a) Persistent current, j, as a function of single im- 
purity strength, £i mp , obtained from the CLDFT for L = 53, 



TV = 26, U/t = 4 and $ 



In (b) we show a typical site 



density profile for a positive single impurity site potential. 



IV. CONCLUSION 

In this work we have presented an extension of the 
BALDA for the one-dimensional Hubbard problem on a 
ring to CLDFT. We have then investigated the response 
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of interacting rings to an external flux both in the ho- 
mogeneous and inhomogeneous case, and we have com- 
pared our results with those obtained by numerically ex- 
act techniques. Our analysis has been confined to the 
metallic limit, i.e. away from half-filling, where the Hub- 
bard model has a metal to insulator transition. In general 
we have found that CLDFT performs rather well in calcu- 
lating both the persistent currents and the Drude weights 
in the homogeneous case. Furthermore a similar level of 
accuracy is transferred to the impurity problem. With 
these results in hands we propose to use CLDFT in the 
study of AB rings where the combined effect of electron- 
electron interaction and disorder can be addressed for 
large rings, so that a numerical evaluation of the vari- 
ous scaling laws proposed in the past can be accurately 
carried out. 



Here D c (n) is the charge stiffness (Drude weight) defined 
as 



_ L d 2 £{n^) 

c 2 



(34) 



In physical terms the Drude weight D c is the real part 
of the optical conductivity a\(w) in the long wavelength 
limit 1611 . 



ax{w) = 2nD c S(w) + a[ es (w) 



(35) 



where we took h = e = c = 1. If we de- 
note £ BA (nBA, ^ba) an d £°(n ,$o) respectively as the 
ground state energy for the interacting system [first term 
in equation ( 18 1] and for the non-interacting one [sec- 
ond term in equation ( 18 )], away from half-filling we will 
write 
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VI. APPENDIX: LOCAL DENSITY 
APPROXIMATION FOR THE CLDFT 

We use the BA solution for the homogeneous part of 
the [equation ^} to estimate the XC energy. Then 
the local approximation is taken, 



i 



[ni,3i\ 



(30) 



Here e xc (= 



£ xc ["j] - 



is the XC energy per site for the ho- 



mogeneous system, which is provided in equation ( 18 1 



The first term of the equation (18) can be calculated ex- 
actly using the BA procedures |59| to obtain the ground 
state energy as a function of n and $. Then the phase 
variable <f> can be eliminated from the ground state en- 
ergy to contain the current via 



3 



(31) 



The complete flux dependence of the ground state en- 
ergy for the Mott insulator phase (n = 1) in the thermo- 
dynamic limit has been shown to be |61j 



£(n, $) - £ (n, 0) = 2Dc ^ (I - cos $) , 
while away from half filling and L — > oo this is 



£ (n, $) - £ (n, 0) 



DJn) 



L 



-<f> 



(32) 



(33) 



£ BA (n BA ,$BA) =£ UA (n BA ,0) + 



D* A (n BA ) 



L 



$ 2 



BA > 



£°(no,$o) =£°(n ,0) + 



D° c (n ) 



L 



and 



/ BA (n BA ,$B A ) =2^^2$ 



o • 



BA 



(36) 



(37) 



The fundamental requirement of the KS mapping is that 
fiBA = no = n and j BA = j° = j while we note that 
$ba = $ and <f> = & s in equation ([18]). By substituting 
equation (36) and the expressions for <f> s and $ obtained 



from equation (|37|) into equation (|I8| one obtains 
£ xc (n,j) 
where 



£ BA (n, 0) - £°(n, 0) - £ H (n) + -A xc (n)j 2 , 

(38) 



A xc (n) 



f 



1 



D°(n) £>ba(„) 



(39) 



Here D®(n) is the non-interacting charge stiffness defined 
as 



„n, s 2i /mr\ 
Dg(n) = - B m( T ) 



(40) 



for L -> oo. Df A (n) can then be obtained in the ther- 
modynamic limit by using [68) 



1 

2^ 



(41) 



where £ c is an element of the dressed charge matrix, 
which is used to describe the scattering between the 
quasi-particles and v c is velocity of the charge excitation. 
Therefore one finally obtains 



e xc (n,j) 



e xc (n,0) + ^A xc (n)j 2 , (42) 
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so that 



^BALDA iP'l i 3l ) = 



<9e xc (n, j) 



dn 



(43) 



n->-ni,3->-3i 



and 



n->ni,j->-ji 



(44) 



J.H. Davies, The Physics of Low-Dimensional Semicon- 
ductor, (Cambridge University Press, Cambridge, 1998). [28 
See for example T. Giamarchi, Quantum Physics in 
One Dimension, First Edition, (Oxford University Press, [29 
2004). 

E. Zipper, M. Kurpas, M. Szelag, J. Dajka and M. Szopa, [30" 
Phys. Rev B 74, 125426 (2006). 

P. Foldi, O. Kalman, M.G. Benedict and F.M. Peeters, [31 
Phys. Rev. B 73, 155325 (2006). 

A. Lorke, R. J. Luyken, A.O. Govorov, J. P. Kotthaus, [32 
J.M. Garcia and P.M. Petroff, Phys. Rev. Lett. 84, 2223 
(1999) [33; 
A. Fuhrer, S. Lusher, T. Ihn, T. Heinzel, K. Ensslin, 
W. Wegscheider and M. Bichler, Nature (London) 413, [34 
822 (2001) 

J.M. Garcia, G. Medeiros-Ribeiro, K. Schmidt, T. Ngo, [35 
J.L. Feng, A.Lorke, J. Kotthaus and P.M. Petroff, Appl. [36" 
Phys. Lett. 71, 2014 (1997). 

A. Lorke and R. J. Luyken, Physica B 256, 424 (1998). [37" 
Y. Aharonov and D. Bohm, Phys. Rev. 115, 485, (1959). 
L.P. Levy, G. Dolan, J. Dunsmuir and H. Bouchiat, Phys. [38 
Rev. Lett. 64, 2074 (1990). 

M. Biittiker, Y. Imry and R. Landauer, Phys. Lett. A 
96, 365 (1983). 

F. Bloch, Phys. Rev. 137, A787 (1965). 

F. Bloch, Phys. Rev. 166, 415 (1968). 

M. Biittiker, Phys. Rev. B 32, 1846(R) (1985). [39 

G. Queeroz-Pallegrino, J. Phys.: Condens. Matter 13 
8121 (2001). 

S. Kirchner, H.G. Evertz and W. Hanke, PRB 59 1825 [40 
(1999). 

R.A. Molina, D. Weinmann, R.A. Jalabert, G. Ingold and [41 
J. Pichard, Phys. Rev. B 67 235306 (2003). 
S.K. Maiti, J. Chowdhury and S.N. Karmakar Phys. Lett. [42 
A 332 497 (2004). 

T. Giamarchi and B. Sriram Shastry, Phys. Rev. B. 51, [43 
10915 (1995). 

E.B. Kolomeisky and J. P. Straley, Rev. Mod. Phys. 68, [44 
175 (1996). 

S. Viefers, P. Koskinen, P. Singha Deo and M. Manninen, 
Physica E 21, 1 (2004). [45 
J. Hubbard, Proc. Roy. Soc. A 276, 238 (1963). 
J. Hubbard, Proc. Roy. Soc. A 277, 237 (1964). 

G. Bouzerar, D. Poiblanc and G. Montambaux, Phys. [46 
Rev. B 49, 8258 (1994). 

H. A. Bethe, Z. Phys. 71, 205 (1931). [47 
E.H. Lieb and F.Y. Wu, Phys. Rev. Lett. 20, 1445 (1968), 

E. H. Lieb and F. Y. Wu, Physica A 321, 1 (2003). [48 
[27] V. Meden and U. Schollwock, Phys. Rev. B 67, 035106 [49 



(2003). 

F.C. Dias, I.R. Pimentel and M. Henkel, Phys. Rev. B 
73, 075109 (2006). 

A. Cohen, K. Richter and R. Berkovits, Phys. Rev. B 57, 
6223 (1998). 

A. O. Gogolin and N.V. Prokof'ev, Phys. Rev. B 50, 4921 
(1994). 

S. Jaimungal, M.H.S. Amin and G. Rose, Int. J. Mod. 
Phys. B 13, 3171 (1999) 

S. Viefers, P. Singha Deo, S.M. Reimann, M. Manninen 
and M. Koskinen, Phys. Rev. B 62, 10668 (2000). 

F. Pederiva, A. Emperador and E. Lipparini, Phys. Rev. 
B 66, 165314 (2002). 

P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 
(1964). 

W. Kohn and L.J. Sham, Phys. Rev. 140, A1133 (1965). 
O. Gunnarsson and K. Schonhammer, Phys. Rev. Lett. 
56, 1968 (1986). 

K. Schonhammer, O. Gunnarsson and R.M. Noack, Phys. 
Rev. B. 52, 2504 (1995). 

K. Capelle, N.A. Lima, M.F. Silva and L.N. Oliviera, in 
The Fundamentals of Electron Density, Density Matrix 
and Density Functional theory in Atoms, Molecules and 
Solids, Kluwer series, "Progress in Theoretical Chemistry 
and Physics," edited by N. I. Gidopoulos and S. Wilson 
(Kluwer, Dordrecht, 2003). 

Density functionals: Theory and Applications, edited by 
D. Joulbert, Springer Lecture Notes in Physics Vol. 500 
(Springer, Berlin, 1998). 

N.A. Lima, L.N. Oliviera and K. Capelle, Euro. Phys. 
Lett. 60, 601 (2002). 

N.A. Lima, L.N. Oliviera and K. Capelle, Phys. Rev. 
Lett. 90, 146402 (2003). 

M.F. Silva, N.A. Lima, A.L. Malvezzi and K. Capelle, 
Phys. Rev. B. 71, 125130 (2005). 

V. L. Campo, Jr. and K. Capelle, Phys. Rev. A. 72, 
061602(R) (2005). 

G. Xianlong, M. Polini, M.P Tosi, V.L. Campo, 
K. Capelle and M. Rigol, Phys. Rev. B. 73, 165120 

(2006) . 

G. Xianlong, M. Rizzi, M. Polini, R. Fazio, M.P. Tosi, 
V.L. Campo and K. Capelle, Phys. Rev. Lett. 98, 030404 

(2007) . 

W. Li, G. Xianlong, C. Kollath and M. Polini, Phys. Rev. 

B. 78, 195109 (2008). 

A. Akande and S. Sanvito, Phys. Rev. B 82 245114 
(2010). 

C. Verdozzi, Phys. Rev. Lett. 101, 166401 (2008). 

S. Kurth, G. Stefanucci, E. Khosravi, C. Verdozzi and 



11 



E.K.U. Gross, Phys. Rev. Lett. 104, 236801 (2010). 
[50] S. Kurth and G. Stefanucci, rarXiv:1012. 4296 
[51] D. Karlson, A. Privitera and C. Verdozzi, 

|arXiv:1004.2264] (2010). 
[52] M. Dzierzawa, U. Eckern, S. Schenk, and P. Schwab, 

Phys. Status Solidi B 246, 941 (2009). 
[53] E. Runge and E.K.U. Gross, Phys. Rev. Lett. 52, 997 

(1984). 

[54] G. F. Giuliani and G. Vignale, Quantum Theory of 
the Electron Liquid, (Cambridge University Press, Cam- 
bridge (UK), 2005). 

[55] G. Vignale and M. Rasolt, Phys. Rev. Lett. 59, 2360 

(1987) . 

[56] G. Vignale and M. Rasolt, Phys. Rev. B. 37, 10685 

(1988) . 

[57] R. Peierls, Z. Phys. 80, 763 (1933). 
[58] M. Graf and P. Vogl, Phys. Rev. B 51, 4940 (1995). 
[59] B.S. Shastry and B. Sunderland, Phys. Rev. Lett. 65, 
243 (1990). 



[60] F. Nakano, J. Phys. A: Math. Gen. 33, 5429 (2000). 
[61] C.A. Stafford and A.J. Millis, Phys. Rev. B 48, 1409 
(1993). 

[62] R.A. Romer and A. Punnoose, Phys. Rev. B 52, 14809 
(1995). 

[63] B.-B. Wei, S.-J. Gu and H.-Q. Lin, J. Phys.: Condens. 

Matter 20, 395209 (2008). 
[64] S. Schenk, M. Dzierzawa, P. Schwab and U. Eckern, Phys. 

Rev. B. 78, 165102 (2008). 
[65] P. Schmitteckert, T. Schulze, C. Schuster, P. Schwab and 

U. Eckern, Phys. Rev. Lett. 80, 560 (1998). 
[66] P. Koskinen and M. Manninen, Phys. Rev. B. 68, 195304 

(2003). 

[67] H-F. Cheung, Y. Gefen, E.K. Riedel and W.H. Shih, 

Phys. Rev. B 37, 6050 (1988). 
[68] N. Kawakami and S.-K. Yang, Phys. Rev. B 44, 7844 

(1991). 



